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Abstract: The production rate of the CP-odd Higgs boson in the Minimal Supersymmet- 
ric Standard Model is evaluated through next-to-leading order in the strong coupling con- 
stant. The divergent integrals are regulated using Dimensional Reduction, with a straight- 
forward implementation of 75. The result is confirmed within Dimensional Regularization 
where 75 is implemented according to the Standard Model calculation of Chetyrkin et 
al. [16]. The well-known Standard Model result is recovered if the masses of the supersym- 
metric particles tend to infinity. 
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1. Introduction 

The Minimal Supersymmetric Standard Model (MSSM) predicts a fundamental CP-odd 
scalar particle A, commonly referred to as the pseudo-scalar Higgs boson. One of the 
most important properties that distinguishes it from its CP-even analogues h and H is 
the absence of tree-level couplings to the electro-weak gauge bosons W and Z. Decay 
and production processes through these particles, which have been shown to be extremely 
helpful for CP-even Higgs searches and studies, are thus very much suppressed. This leaves 
associated ttA and bbA production as well as the loop-induced gluon fusion process as the 
most important production modes of a pseudo-scalar Higgs boson at the Large Hadron 
Collider (LHC) (for a recent review, see Ref. [1]). In this paper, we present the evaluation 
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of the inclusive gluon fusion cross section through next-to-leading order (NLO) in the strong 
coupling constant. 

Despite the fact that the tree-level gg4> coupling vanishes {(j) £ {h, H, A}), gluon fusion 
in general has a comparatively large cross section because of the high gluon luminosity 
at the LHC, and the large top- Yukawa coupling. In the limit where squarks and gluinos 
are decoupled, QCD corrections have been evaluated through NNLO [2-6] by employing 



an effective Lagrangian for heavy top quarks (cf. Sect. 3^); they have been shown to be 
numerically significant but perturbatively well-behaved. 

A precise determination of the gluon fusion cross section at the LHC could yield sensitivity 
to as yet undiscovered particles that may mediate the gg4> coupling apart from the top 
and bottom quarks. Within the MSSM, for example, top squarks can play a significant role 
if they are lighter than around 400 GeV (recall that the Yukawa coupling for squarks is 
typically proportional to rrig rather than m?). In the case of CP-even Higgs bosons, such 
effects occur already at leading order (LO) (i.e., 1-loop). The NLO result was obtained in 
Ref. [7] within an effective theory for top, stop, and gluino masses much larger than the 
Higgs mass. 

For the CP-odd Higgs boson, however, squarks do not affect the ggA vertex at 1-loop level 
due to the structure of the Aqq coupling, as will be shown below. The two-loop effects are 
thus expected to have a larger influence than for the CP-even Higgs production. What adds 
to this is that, in the limit of large rrig, the quark mediated contribution to the ggA vertex 
does not receive any QCD corrections owing to the Adler-Bardeen theorem. As we will 
show, the only 2-loop QCD effects to this coupling are due to mixed gluino-quark-squark 
diagrams in this limit, leading to a potentially increased sensitivity to the gluino mass. 
Of course, the heavy quark limit is not applicable for the bottom mediated gluon-Higgs 
coupling which does receive QCD corrections [8-10]. 

The calculation involves a technical issue that deserves special mention, namely the imple- 
mentation of 75 within Dimensional Reduction (DRED). Its mathematically consistent and 
practically feasible formulation has been a subject of interest for many years now [11] (for 
recent developments concerning DRED, see Ref. [12,13]). Here we adopt an approach close 
to the prescription of Refs. [14, 15]. In addition, we perform the calculation using Dimen- 
sional Regularization (DREG) with 75 implemented as in Ref. [16] and find full agreement 
with the result obtained in DRED. 

Another important consistency check for our calculation is obtained from the SM limit: even 
for an infinitely heavy supersymmetric (SUSY) spectrum, the SUSY diagrams give a non- 
vanishing contribution which exactly cancels the top mass counter-term contribution arising 
from the LO SM diagram. The well-known result of vanishing higher order corrections to the 
quark-mediated ggA coupling, as required by the Adler-Bardeen theorem, is thus recovered 
in a non-trivial way. Combining this observation with the considerations of Ref. [13], it 
seems possible that in a supersymmetric theory the formal difficulties of DRED do not pose 
serious technical problems in practical calculations. 
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The outhne of this paper is as follows: In Sect. ^ we introduce our notation and quote the 
LO result for the partonic process gg — A. In Sect.^, the effective Lagrangian underlying 
our calculation is introduced and the treatment of 75 is discussed. Sect.|^ outlines the 
method of the calculation and discusses the general structure of the result. It also provides 
analytic formulae in some limiting cases. The numerical influence of the NLO terms is 
discussed in Sect.|5|. In Sect.|6| our findings are summarized and an outlook on possible 
extensions of this work is given. 



2. Notation and Leading Order Result 
2.1 Lagrangian 

We write the underlying Lagrangian in the following form: 

^ = -C-QCD + -C-qA + ^SQCD + ^qA , (2-1) 

where 

2 2 

Aa = ^5Z^5^™^, CqA = Y.Y.'^~9tijmjA. (2.2) 

q q i=l 

The sums run over all quark flavors, is the mass of the quark q and v ~ 246 GeV the 
Higgs vacuum expectation value; the coupling constants g^ and g^^^ will be defined below. 
^QCD denotes the full QCD Lagrangian with six quark flavors, while >Cqcd + ^SQCD is 
the supersymmetric extension of £qcd within the MSSM, i.e., >CsQCD incorporates kinetic, 
mass, mixing, and interaction terms of all the squarks and gluinos. Since we will be 
concerned with higher orders in the strong coupling only, A does not appear as a 
dynamical field and does not require a kinetic or a mass term. 

qi,q2 denote the squark mass eigenstates which are related to the chiral eigenstates qL,QR 
(the superpartners of the left- and the right-handed quark q) through 




cos ( 



sm ti„ cos 




(2.3) 



yq i^uo wq 

The coupling constants relevant for the following discussion are^ 
g^ = tan /3 , gf = cot /? , gA^ = = ^ 



m 



2 ^2 



~9^i2 = -~9t2i = \ 2 sin 2dt cot 13 + (cot^ /? + 1) 



The effect of quarks other than bottom and top can be neglected because of their small 
Yukawa couplings. Bottom and sbottom effects are typically of order (m^/M^) tan^ /? 



^For a more complete collection of Feynman rules, see Ref. [17], for example. 
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relative to the top and stop effects. In this first analysis of NLO squark effects to pseudo- 
scalar Higgs production, we neglect the effect of sbottom quarks; they require a further 
extension of our approach. Our results are thus strictly valid only for not-too-large values 
of tan/3. Pure bottom effects, on the other hand, will be included through NLO in our 
numerical analysis. 



2.2 Hadronic cross section 

The cross section for the hadronic process pp ^ A + X at a center-of-mass energy -^i is 
determined by the formula'^ 

'^{z)= ^ dxi dx2 ipi{xi)ipj {x2) aij i j, z := — , (2.5) 

where '^i{x) is the density of parton i inside the proton, dij is the cross section for the 
process ij ^ A + X, where, as indicated in Eq. ( |2.5D , i and j are parton labels. For our 
numerical analysis, we will use the MRST parton density sets throughout this paper [18,19]. 

Sample diagrams for ij=gg, gq,qq that contribute to the inclusive Higgs production rate 
at LO and NLO are shown in Fig. ||. 



2.3 Standard Model limit and leading order result 

In the limit where all SUSY masses tend to infinity (denoted in the following by MsusY 
oo), the Lagrangian of Eq. ( p.lD reduces to 

C^^ = Cqcd + ■ (2.6) 



This will be called the "Standard Model (SM)" limit in what follows, despite the fact that 
the SM does not contain a CP-odd Higgs boson. Consequently, diagrams without any 
squark and gluino lines will be called "SM diagrams" (e.g. Fig.|l|(a)-(f)), while "SUSY 
diagrams" contain at least one propag ator of a SUSY particle (Fig.^ (g)-(j)). 

Note that due to the antisymmetric structure of the g^^j , see Eq. (|2.4D , there are no SUSY 
diagrams at the one-loop level. The LO result for the partonic process gg ^ A is thus 
determined solely by diagrams of the type shown in Fig. |l](a). It reads explicitely [20-23]: 



6{1 - x) + 0{al) 



256^2 V vr 



q£{t,h} 



(2.7) 



= 9n rgf{Tg) , Tg 



Arrig 
Ml 



^The modifications for pp collisions are obvious and shall not be pointed out explicitely in this paper. 
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(c) 



(d) 




.-A 





(g) 



■A 

(h) 



(i) 




(j) 



Figure 1: Diagrams contributing to the process pp A + X in the MSSM at (a) LO and (b)~(j) 
NLO. According to the definition of the main text, (a)-(f) are referred to as "SM diagrams", (g)-(j) 
as "SUSY diagrams" . The straight-curly hne connecting quarks and squarks represents a gluino 
(denoted by g in the text). 



where 



fir) 



arcsm 



In 



1-VT 



ITT 



T> 1. 



r < 1. 



(2.8) 



Throughout this paper, ag denotes the strong coupHng constant, renormaUzed in the MS 
scheme for QCD with five massless quark flavors. 

3. Higher orders 

3.1 Standard Model result — Effective Lagrangian 



The SM contributions, based on £sm of Eq. ( p.q) , are known through NLO in terms of 
1-dimensional integral representations [8,9], implemented in the program HIGLU [24]. They 
include both top and bottom quark loops as shown in Fig. ||(a)-(f) {q G {b,t}). 

It has been shown in Ref. [9, 10] that the NLO top quark contributions are well approxi- 
mated by the formula 

,NLO 



cr, 



NLO 



LO 



where 



LO 



(3.1) 
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(Too is the cross section evaluated in an effective theory, obtained by integrating out the top 
quark: 

A 



mt'^MA £SM 



V L 



■SM , 



The operators are defined as 



+ C 



(5) 
QCD 



(3.2) 



(3.3) 



Here, {q} := {d,u, s,c,b} denotes the set of hght (in our case massless) quark fields and 

(5) 

GJ^^ the gluon field strength tensor, /^qcd i^ Standard QCD Lagrangian with five quark 
fiavors. 

The Wilson coefficients Cf^ and Cf^ have been calculated through NNLO in Ref. [16]. 
Cf^ contributes only at NNLO and shall not be discussed any further in this paper. The 
resuh for Cf^ IS 



. cot /3^ + 0{at) 



(3.4) 



as it had been previously suggested in Refs. [8, 25, 26] on the basis of the Adler-Bardeen 
theorem [27]. 

Oi generates vertices which couple two and three gluons to the pseudo-scalar Higgs boson 
(the ggggA-veitex vanishes according to the Jacobi identity of the structure functions of 
SU(3)). Sample diagrams that contribute to the NLO cross section in the effective theory 
of Eq. ( |3.2| ) are shown in Fig. ||. The full set has been calculated through NLO in Ref. [8,9], 
and through NNLO in Ref. [3,4,6]. 




(a) 








(d) (e) (f) 

Figure 2: Diagrams contributing to the gluon fusion process in the effective theory at NLO. 



3.2 SUSY contributions 

We construct again an effective theory, this time by integrating out not only the top quark, 
but also all the SUSY particles. Since the remaining degrees of freedom are the same as in 
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the Standard Model case of Sect. 3.1, the effective Lagrangian has exactly the same form 



as in Eq. (|]|), only with Cf^ and Cf^ replaced by Ci and C2. 

As in the SM case, C2 is zero through NLO. It remains to determine Ci within the framework 
of Eq. ( |2.1| ) through NLO in a<j. The calculation of the SM case in Ref. [16] was done 
using Dimensional Regularization as it is most convenient for loop calculations in QCD. 
It is known that DREG breaks SUSY and requires SUSY-restoring finite counter terms in 
general. In the current case, however, these are not required due to the absence of SUSY 
contributions at LO. 

On the other hand, one may apply Dimensional Reduction [28] in order to regularize the 
divergent integrals. This is done by setting D = 4: after contracting all (Lorentz and 
spinor) indices [29,30], while the loop integrals are subsequently evaluated in D = 4 — 2e 
dimensions. 

A particularly subtle issue is the treatment of 75 which occurs in the Higgs-quark-anti- 
quark and the gluino-quark-squark vertex. We anticommute these 75 matrices and use 
7I = 1 until only one of them remains in the Fermion trace. Working in DREG, we 
may follow the calculation of Ref. [16] which adopted the method of Ref. [14,31,32] for 
the treatment of 75. This involves a finite counter term in order to restore gauge 
invariance [33]. Since there are no SUSY diagrams at LO, is given by the well-known 
SM expression quoted in Ref. [16]. 

For the DRED calculation, on the other hand, we simply set [14,16,32] 

75 = |yW7'^7V7'^. (3.5) 

It has been shown that the occurrence of the Levi-Civita symbol e^upa can lead to incon- 
sistencies when implemented in DRED [11]. We circumvent them by keeping the genuinely 
4-dimensional Levi-Civita symbol Ef^upa uncontracted until after the renormalization pro- 
cedure, in close analogy to Refs. [16,32]. As opposed to the calculation in DREG [16], 
however, working in DRED by definition does not involve terms of 0{e = 2 — D/2) that 
require the above-mentioned finite counter terms in order to restore gauge invariance. We 
have verified that the scalar and pseudo-scalar one-loop quark vertices are identical in the 
chiral limit. Thus, the formulae of Ref. [16] for projecting onto the coefficient function Ci 
can be translated to the DRED case by setting D = 4 and ignoring the finite renormalization 
of the pseudo-scalar current (i.e., Z^ = 1). 

We will explicitely demonstrate that the SM result is recovered (through 0{ag)) by making 
all the superpartner masses infinitely heavy. It turns out that the SUSY diagrams lead to 
terms that do not vanish as MgusY 00. However, they cancel with the counter terms of 
the SM diagrams to give the correct SM result. 

The generalization of this approach for the implementation of 75 within DRED through 
higher orders is very desirable, of course, because it greatly simplifies precision calculations 
in SUSY models as they might be required by future experimental data. 
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4. Results 



The evaluation of the diagrams proceeds in complete analogy to Ref. [7]. At NLO, one needs 
to evaluate massive 2-loop diagrams with vanishing external momenta. This is possible in 
a fully analytic way with the help of the algorithm of Ref. [34]'^. As a cross check, we 
also calculated the diagrams in the limit rrit ^ m^^ <^ rrif^ <^ rrig by using automated 
asymptotic expansions [36] and found agreement with the corresponding expansion of the 
analytical result. 

Subsequently, the bare coupling constant within the SUSY-QCD theory is transformed to 
its renormalized 5-flavor QCD expression in the MS scheme as described in Ref. [7]. 

In DRED, and using the prescription for the treatment of 75 described above, we find the 
following contribution of all SM diagrams to the coefficient function: 

The second term in the curly brackets arises from the MSSM expression for the quark 
mass counter term. Using the SM expression instead, we recover the well-known result of 
Eq. dU). The SUSY diagrams add up to 

^lIsUSY-dias = "I^ ^{ 3 V| + ^^jh^ + ' ^ ^ ^"^^"'"^^^ ' ^'-'^ 

such that indeed the SM limit is recovered as the SUSY masses are decoupled: 

=(SM 



-«,r-col'^'lsM-=.,» + C^l|sUSV-.,J' (4.3) 

The general NLO result for Ci can then be cast into the following form: 

Ci = --^cot/?|l + — ^^^^(cot/? + tan/3) f{mt,mr,mr,mg)\+0{al). (4.4) 
iovr 1^ vr rrig ^ J 

A few observations may be worth pointing out: 



Eq. ( |4.4D immediately shows that the NLO corrections to the coefficient function will 
get more important w.r.t. the LO term for large values of tan/3 and //suSY- On the 
other hand, large values of tan/3 increase the importance of bottom contributions 
even more, so that they usually obscure the squark effects. 

The 0{a1) corrections are proportional to ^susYi while the squark mixing angle 9t 
(or, equivalently, the trilinear coupling At) drops out in the final result. This is due 
to the axial U(l) Peccei-Quinn symmetry of the SUSY potential [23,37] whose explicit 
breaking by the /xsusY term violates the Adler-Bardeen theorem. It provides another 
consistency check of our calculation. 



^We are indebted to M. Steinhauser for providing us with his implementation of this algorithm in the 
framework of MAT AD [35] . 
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• As a result of the absence of mass terms at LO, the expHcit form of / is independent 
of the mass renormahzation scheme. 

• Eqs. (^]^) and ( |4.2| ) are unchanged in DREG if the SM diagrams are multiphed by 
the finite renormahzation constant (see Sect. and Ref. [16]). 

The general result for f {mt, m^_^, m^^, nig) is too voluminous to be displayed here. Instead, 
we implemented Ci in the program evalcsusy.f [7]^ and analyze its behavior numerically 
as shown in Sect.|5|. 

For practical purposes, it might be useful to provide the analytical expression of the function 
/ for some limiting cases. Defining two mass scales m <C M, and 

X = ^ , l. = ln(x) , (4.5) 



we find 



5 9 

f{m,m,m,m) = —-- 52 = -0.894176 ... , (4.6) 

18 2 

f{m,M,M,M) = ---x(—-—l^]-x^(^ + — l^]+..., (4.7) 
•'^ ' ' ' ^ 3 V54 36 y V1800 60 7 ' v / 

22 /71 13 1\ 

f{m,m,m,M) = - + a; f- -- C2 + y g ^'j +... , (4.8) 

3 3, /II „ . 29 , 3 2 



f{M,m,m,m) = - + - + 2; (^y + 3 C2 + y + - ^^J + • • • , (4.9) 

13 /22 31 3 \ 

fim,M,M,m) = - _ x - (^y + 3 C2 + y + - j + • • • , (4.10) 

2 /I8I 13 \ 

f(m,m,M,M) = 2; + ]+..., (4.11) 

where the dots denote higher orders in x, and 



Insertion of f(m,m,m,M) into Eq. (4^) shows that, in contrast to scalar Higgs produc- 
tion [7,38,39], the result for Ci is well-behaved as rUg — > 00 at finite m. This is because 
there are no squark contributions at LO which could affect the renormahzation of the NLO 
terms. The same holds for — > 00 as can be seen from the form of f{m, M, M, m). Note, 
however, that the result is logarithmically divergent in the case where mt — > 00 with the 
other masses fixed, see f{M,m,m,m). This is due to the linear mass dependence of the 
top Yukawa coupling and can be observed already in the pure SM contributions at higher 
orders in a,. 



' The program is available from http : //www-ttp .physik .uni-karlsruhe . de/Progdata/ttp04/ttp04-19 . 



-9- 



5. Discussion 



5.1 Coefficient function Ci 



We will now discuss the numerical effect of the newly evaluated terms. First, we investigate 
the dependence of Ci on the squark and gluino masses. To this aim, it is convenient to 
consider the leading and the NLO term separately: 



According to Eq. 



(5.1) 



cf^ = cot 13. 



(5.2) 



Fig. ^ shows ^ as a function of rrig for various values of m^^ and m^^. In Fig.|^, on the 
other hand, we consider it as a function of one of the stop masses {mi^) while fixing rrig and 
the other stop mass {mi^) at a few representative values; A and B assume the values 1 and 
2, depending on whether m^^ is larger or smaller than m^^ . Note that is symmetric in 
rrij^ and m^^ since it does not depend on the squark mixing angle 9t as mentioned above. 



Finally, in Fig. ^, 



m 



= mi is varied for certain choices of rrig. 



The general structure is quite similar in all figures: c\ is of the order of —0.5 for moderate 
values of the SUSY masses. As they increase, c^l^ tends to zero in agreement with the 
general discussion above. In order to demonstrate this behavior more clearly. Fig. ^ adopts 
the same parameters as Fig. H, but extends up to very large values of the gluino mass. 
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Figure 3: c"^^ as a function of the gluino mass for mt = 173 GeV, /isusY = 150 GeV, tan/? = 3, 
and (TOt^,TOtJ ^ (200, 200)7(200, 400)7(200, 600)7(400, 600) GeV [long7. . . 7short dashes]. 
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Figure 4: c'^^ as a function of m^^ for mt — 173 GeV, /isusY = 150 GcV, tan /3 = 3, = 400 GcV, 
and mi^ = 200/300/400/600GeV [long/. . . /short dashes]; [(A, B) = (1, 2) if m^^ < TOj^ , otherwise 
(AS) = (2,1)]. 
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Figure 5: c^^-* as a function of TOj 
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mi for mt = 173 GeV, /Ususy = 150 GeV, tan/3 = 3, 



200/400/600/800 GeV [long/. . . /short dashes]. 



For the discussion of the numerical effects on the cross section in the following section, we 
adopt a scenario similar to the m™^ benchmark scenario [40] , corresponding to 

ruf^ = 826 GeV , mj^ = 1172 GeV , m-g = 800 GeV . (5.3) 

In the original definition of this scenario, fisvSY was set to —200 GeV; here, however, we 
will allow for a variation of this parameter between ztlTeV. In addition, we will assume 
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Figure 6: Same as Fig.0, but extended abscissa in order to demonstrate tlie asymptotic behavior 
or ' . 

tan/3 = 3; the top mass is set to rat = 173 GeV.^ 

The 2-loop term of the coefficient function is then given by 



-0.392 



^SUSY 

ITeV 



(5.4) 



5.2 Effects on the cross section 



The inclusive NLO hadronic cross section for pseudo-scalar Higgs production receives con- 
tributions from the subprocesses gg — > ^(-l-^), qg — > Aq^ and qq — > Ag. At this order, 
squarks only affect the process gg ^ A, because of the antisymmetric structure of the 
qtqjA coupling, cf. Eqs. (^.21) , (|2.4D . Recalling the notation of Eq. (|2.5| ), we write 

^ggix) = atbix) + Aa^^ix) + Aaif^ix) + 0{at) , (5.5) 

where atb denotes the contributions arising from the top and bottom mediated gluon-Higgs 
couplings. It can be evaluated through NLO for arbitrary top, bottom, and Higgs masses 
with the help of the FORTRAN program HIGLU. The same is true for the qg and the qq 
sub-processes (see Figs.|l](e) and (f)). 

A(T^j and Act^^ are the effects arising from the interference of stop-induced with top- and 
bottom-induced amplitudes: 

Aa^g ~ Re (^Mf^M^^')) , q G {b,t} . (5.6) 

^ We remark that, for Ma ^ 200 GeV, this choice of parameters leads to a conflict between the theoretical 
value of the light Higgs mass (evaluated with FeynHiggs [41] , for example) and its current experimental 
lower limit; nevertheless, for the sake of generality, we will vary Ma between 100 and 300 GeV in our 
numerical analyses below. 
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Ai'^^ is expressed in terms of Aq from Eq. (2.7), while for we use the expression 

evaluated in the effective theory to obtain 



TT 



^^k(^) = f ^) R« (~^? Mr,) ) S{1 - x) . (5.7) 

Note that -M^^^ has a branch cut at Ma = 2mt. Thus, we expect our result to be valid 
for Ma < '^mt- Recall, however, that in the SM case, the heavy top limit still provides an 
excellent approximation for Higgs masses much larger than 2mt [9,10,42,43]. 

To study the numerical effects, we consider the modified m™^^ scenario quoted in Eq. (^). 
Fig. shows the relative size of the top squark effects to the total NLO cross section. We 
note that even for |/isusY| = 1 TeV, they hardly exceed 4%. 

Fig. ^ shows separately the effects of the top-stop and the bottom-stop interference terms, 
again relative to the total NLO cross section, for ^ = 1 TeV. They are of similar order in 
magnitude, but opposite in sign, thus cancelling each other to a certain extent. 

Finally, Fig. |9| shows the inclusive cross section through NLO for the modified m™''^ scenario 
defined in and below Eq. (|5.3|) , including effects of top and bottom quarks (solid line) as 
well as top squarks (dashed lines) for //susY = ±lTeV. 
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Figure 7: Effects of the top-stop and bottom-stop interference terms relative to the total NLO 
cross section for the scenario defined in and below Eq. ( ^.3[ ). The numbers above the graphs denote 
the value of /isusY in GeV. 
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Figure 8: Relative effect of the top-stop (dotted) and bottom-stop (dashed) interference terms, as 







well as their sum (solid), for the scenario defined in and below Eq. (5.3) 



6. Conclusions 



The corrections to the effective Higgs-gluon coupling for a pseudo-scalar Higgs boson have 
been evaluated in the MSSM through first order in the strong coupling constant ag, taking 
into account effects of top and bottom quarks as well as top squarks. The NLO corrections 
are proportional to /isusY as expected from the symmetries of the SUSY potential. The 
numerical effects were studied within a specific SUSY scenario, derived from the m™^^ 
scenario of Ref. [40]; further studies can be performed easily using the publicly available 
numerical routine evalcsusy .f (see footnote on page P). 

The calculation also addresses a technical issue, since it involves the 75 matrix in a non- 
trivial way. In analogy to Ref. [16], we avoided the contraction of the Levi-Civita symbol 
with D dimensional quantities. We argued that the calculation in DRED does not require 
finite counter terms as opposed to the DREG approach, provided the underlying theory 
is supersymmetric. The result was confirmed by evaluating the diagrams in DREG and 
including the proper counter term for 75. 

It would be interesting to investigate these observations in more detail, in particular to 
prove the validity of our implementation of 75 within DRED in a rigorous way. Further 
corroboration could be obtained from its application at second order as- This corresponds 
to the evaluation of Ci at three loops for which the technical tools are in principle avail- 
able [44]. 

From the phenomenological point of view, inclusion not only of bottom but also sbottom 
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Figure 9: Total cross section at NLO including top and bottom effects (solid). Stop effects (dashed) 
are shown for the scenario defined in and below Eq. ( |3.3[ ) for two choices of hsvsy- 

quarks would be desirable. This could be done in the heavy sbottom limit using asymptotic 
expansions of Feynman diagrams. Moreover, one could calculate the photonic decay rate 
of the CP-odd Higgs boson in a very similar fashion. There, however, one does not need to 
rely strictly on the effective Lagrangian, but could evaluate the first few terms of a Taylor 
expansion in M^/M^ along the lines of Ref. [45], where M G {mt,mj:,mg}. 
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